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1 ABSTRACT 

The mass reduction of satellite solar arrays results in signifi- 
cant panel flexibility, so possibly striking one another dynam- 
ically leading ultimately to structural damage. To prevent this, 
rubber snubbers are mounted at well chosen points of the struc- 
ture and they act as one sided linear spring; as a negative con- 
sequence, the dynamic of these panels becomes nonlinear. The 
finite element approximation is used to solve partial differen- 
tial equations governing the structural dynamic. The models 
are validated and adjusted with experiments done in the ISVR 
laboratory, Southampton university. 




significantly if the imposed displacement is replaced by a pe- 
riodic force excitation, the configuration used is easiest to be 
realized from a technical point of view as the rig is very simple 
to built [see figure[2] 

The present study is to simulate the behavior of a beam which 
strikes a snubber under a periodic excitation. As the interest is 
to deal with the first three eigen frequencies, the beam is mod- 
eled by ten finite elements. 

The results are presented and compared in the frequency do- 
main, the FFT is applied to the predicted and the experimental 
displacement of the last node of the beam finite element. The 
mass effect of the force transducer is taken in account in the fi- 
nite element code. Note that no treatment is done or used from 
the acquisition system as the problem is nonlinear, the transfer 
function and the other functions use linear assumptions. The 
time signal is taken and the processing is done using an exter- 
nal software (Scilab [8]). The numerical predictions compared 
to experimental results show very good agreement. 



Fi guxe 1 Solar array of the satellite under a test on a shaker 




2 INTRODUCTION 

The study of the total dynamic behavior of solar arrays in a 
folded position with snubbers are so complicated, that to sim- 
plify, a solar array is modeled by a clamped-free Bernoulli 
beam with one-sided linear spring. This system is subjected 
to a periodic excitation force. 

The real state of the problem is close to a beam with unilateral 
contact subjected to a periodic imposed displacement of the 
base, but the dynamical behavior of the system do not change 




Figure 2: The rig used for the experiments: a linear clamped-free 
beam in contact with a rubber 
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3 THE MODEL 



u(x,t): verticale displacement 



t-J.\L'iliilio[] Force 



Clamped— Free beam with a unilateral elastic contact 



Figure 3: beam system with an unilateral spring under a period: 
excitation 
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The beam motion with a snubber can be modeled as: 

pSu{x,t) +EIu {iv) (x,t) = F(t)5 Xo - k T u{x x ,t)- (3.1) 
u(0,t) = 0,d x u(0,t) = 

u(x, t) if u < 



u(x, t)_ = 



ifu>0 



(3.2) 



The classical Hermite cubic finite element approximation is 
used, it yields an ordinary differential system in the form: 



Mq + Kq = -[k r (q ni )-]e ni +F(t)e r . 



(3.3) 



Where M and K are respectively the mass and stiffness 
assembled matrices, q is the vector of degrees of freedom of 
the beam, qi — (ui, d x Uj), i = 1, n, where n is the size of 
M, ni and are the indices of the nodes where the spring 
and the excitation force meet the beam respectively. Numer- 
ical time integration was performed using ODE numerical 
integration for 'stiff problems, package ODEPACK is called 
and it uses the BDF method. 

The use of a small electrodynamic shaker yields a technical 
problem due to the reaction of the beam, the input force F(t) 
can not be a simple sine wave. To deal with this problem, a 
force transducer is fixed between the shaker and the beam. 
The numerical codes use the actual force time signal coming 
from the acquisition system which is periodic. Figure [5] 
shows the input force signals with its spectrum contents. 



4 EXPERIMENTAL AND NUMERI- 
CAL RESULTS 

Figures J2] and (3) give a clear idea of the experiment setup. In 
the table, all the parameters used in the simulations and the ex- 
periments are shown as well as the first three predicted eigen 
frequencies of the linear beam without the spring. The inte- 
gration time is fixed at t — Is for all the simulations and the 
experiment sequences which is five times the greater period of 



the system. 

Figure [7] shows the FFT of the numerical and the experi- 
mental displacements for 32HZ and 12AHz, the height of the 
peaks are normalized by the maximum; the predicted frequen- 
cies found are exactly the same measured for a large number of 
harmonics. However, a small shift in the height of these pikes 
appears from the fifth harmonic, this is maybe due to the low 
number of finite elements used to model the beam; some inves- 
tigations are in progress to understand this aspect. Other tests 
with random excitations have shown good agreement. The ef- 
fect of the unilateral contact is underlined, the input frequency 
is split into its all harmonics. 

Figures [5] and [6] show the input excitation force, it is clear 
from the time signal and from its frequency content that this 
force is not a single sine wave, but it is periodic. 
Figure j4] shows the predicted displacement for an excitation 
of 32Hz and \2AHz, the displacement is almost positive, this 
is due to the high stiffness of the spring but the time response 
is still periodic. 

Note that the frequencies 32Hz and \2\Hz are chosen and 
presented herein because they represent the first two nonlinear 
normal modes of the nonlinear normal mode, they are com- 
puted numerically by a numerical sweep test. Non linear nor- 
mal modes NNM of a nonlinear differential systems are an 
extension of the well known linear normal modes, another in- 
terest of the authors with some papers in progress. 



beam length 


beam width 


beam thickness 


beam's Young modulus 


0.35m 


0.0385m 


0.003m 


69 X 10 9 N/m 2 


Spring stiffness 


1 linear eigen frequency 


2 linear eigen frequency 


3 linear eigen frequency 


57.14 KN/m 


19.97Hz 


122.2Hz 


318.8Hz 




Figure 4: The predicted displacements for an excitation of 
32 Hz and 124 Hz, the high stiffness of the rubber yields 
almost a positive displacement 
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Figure 5: The time signal and its FFT of the input force and Figure 6: The time signal and its FFT of the input force and 
the experimental acceleration for an excitation of 32Hz the experimental acceleration for an excitation of 124Hz 



References 



FFT of the numerical and experimental displacement, Excitation =32 Hz. 




Frequency F/Excitation 



FFT of the numerical and the experimental displacement, Excitation =124 Hz. 




Frequency F/Excitaton 



[1] J. H. Bonsel, R. H. B. Fey. and H. Nijmeirjer: Application of a 
Dynamic Vibration Absorber to a Piecewise Linear Beam System. 
Nonlinear Dynamics 37: 227-243, 2004. 

[3] D.Jiang,C.Pierre, S.W. Shaw: Large-amplitude non-linear nor- 
mal modes of piecewise linear systems. Journal of sound and vi- 
bration 272 (2004) 869-891 

[3] Rob H.B. Fey and Berend Winter, Jaap J. Wijker 

Sine sweep and steady - state response of a simplified solar array 
model with nonlinear support. DETC99/VIB-8027, 1999 ASME 
Design Engineering Technical Conferences, Las Vegas, Nevada, 
USA. 

[4] R.H.B. Fey and F.P.H. van Liempt : 

Sine sweep and steady-state response of a simplified solar array 
model with nonlinear elements. International conference on Struc- 
tural Dynamics Modelling, 3-5 june 2002, pp.201-210. 

[5] M Geradin Theorie des vibrations: applications a la dynamique 
des structures. Masson Paris 2eme edition 1996 

[6] Stephane Junca et Bernard Rousselet: Asymptotic Expansions of 
Vibrations with Unilateral Contact, proceeding of the GDR 2501- 
Anglet 2-6 Juin 2008. 

[7] Richard H. Rand:Lecture Notes on Nonlinear Vibrations, Dept. 
Theoretical Applied Mechanics Cornell University Ithaca NY 
14853 

[8] Scilab [www.scilab.org] 

[9] E.L.B. Van De Vorst, M.F.Heertjes D.H. Van Campen and A. De 
Kraker and R.H. B. Fey: Experimental and numerical analysis of 
the steady state begavior of a beam system with impact. Journal of 
Sound and Vibration 1998 212(2), 210-225 



Figure 7: dashed line: FFT and the experimental displace- 
ments, solid line: FFT of the numerical displacements. Two 
different excitation: 32Hz and 124/Jz 
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Numerical displacement at 100 Hz, Tf=1s 
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